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Abstract Quiescent prominences are thin slabs of cold, dense plasma embedded in the 
much hotter and rarer solar corona. Although their global shape is rather irregular, they 
are often characterised by an internal structure consisting of a large number of thin, parallel 
threads piled together. Prominences often display periodic disturbances mostly observed in 
the Doppler displacement of spectral lines and with an amplitude typically of the order of 
or smaller than 2-3 km s , a value which seems to be much smaller than the characteristic 
speeds of the prominence plasma (namely the Alfven and sound velocities). Two particular 
features of these small amplitude prominence oscillations is that they seem to damp in a 
few periods and that they seem not to affect the whole prominence structure. In addition, 
in high spatial resolution observations, in which threads can be discerned, small amplitude 
oscillations appear to be clearly associated to these fine structure constituents. Prominence 
seismology tries to bring together the results from these observations (e.g. periods, wave- 
lengths, damping times) and their theoretical modelling (by means of the magnetohydrody- 
namic theory) to gain insight into physical properties of prominences that cannot be derived 
from direct observation. In this paper we discuss works that have not been described in previ- 
ous reviews, namely the first seismological application to solar prominences and theoretical 
advances on the attenuation of prominence oscillations. 

Keywords waves, MHD, Sun: prominences, Sun: filaments, Sun: atmosphere, Sun: 
Corona, Sun: oscillations 



1 Introduction 

The solar corona is often populated by peculiar objects, dense clouds of cold plasma in- 
explicably floating tens of thousands of kilometres above the photosphere. Such objects 
are routinely seen during solar eclipses, when they can be easily distinguished by their red 
glow, but they can also be unveiled with the help of filters, such as Ha, devised to observe 
the chromosphere (Fig. [TJ. To put it in simple words, these objects (usually called promi- 
nences or filaments) are like chunks of chromospheric gas defying the downward pull of 
gravity and staying in a place higher than the one that apparently corresponds to their large 
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Fig. 1 Left: Ha image of the full solar disk showing the dark, long and thin structures known as solar promi- 
nences or filaments. Right: high-resolution image of the central part of the filament in the white rectangle of 
the left panel. The fine structure of the prominence, in the form of narrow threads, is clearly visible. Images 
taken from |Lin et iu~1l2007) . 



density. This is not the only enigma around solar prominences. In contrast with the MK tem- 
perature of the surrounding corona, prominences remain at a comparatively cool 10,000 K, 
which prompts one to ask what prevents the mechanisms that heat the corona from also 
raising the temperature of prominences and consequently dispersing them. Other pieces of 
the prominence puzzle concern their beginning and end: first, one may wonder not just how 
prominences form but also why they are born in an adverse environment. Second, despite 
their internal dynamics, prominences that have been stable and healthy for weeks suddenly 
die in an spectacular eruption. Hence, like all living beings these creatures are born, grow, 
(do not reproduce) and die, but unlike many living beings the processes shaping the lifetime 
of prominences are largely unknown. Nevertheless, the intervention of a decisive element is 
quite clear: the magnetic field, that is central to all the mentioned processes. 

The reason of our limited insight into the nature of prominences probably has three 
causes (Vial 1998): there is no such thing as a canonical prominence, but a wide range 
of parameters is observed in different objects; no prominence has a uniform structure, but 
they are made of thin threads (or fibrils) and, in addition, different parameter values can be 
detected in different parts of a prominence; and no structure is really isolated, so it is neces- 
sary to understand the physics of the prominence-corona interface, the effect of the coronal 
radiation field (e.g. Anzer and Heinzel 2005 ) and to trace the magnetic fields permeating 
the prominence to their origin at the Sun's surface (e.g. Lin et al.]|2005|l . Our knowledge 
about prominences has been well reviewed by Tandberg-Hanssen (1995); Martin (1998); 
|Patsourakos and Vial| ( |2002"| ), where ample information on the topic can be found. 

One of the first studies about prominence oscillations was carried out by Harvey ( 1969| >, 
who noted that in a sample of 68 non-active region prominences, 31% of the objects pre- 
sented no significant velocity change along the line-of-sight, 28% showed apparently ran- 
dom line-of-sight velocity variations and 41% presented a definite oscillatory behaviour. 
Similar results were obtained for a set of 45 active region prominences. There are several 
reasons that may lead to the absence of periodic variations in some prominences, e.g. the 
velocity amplitude or its projection along the line-of-sight are too small to stand above the 
instrumental noise level; or the prominence material does not actually oscillate at the time 
the observations are performed; or the light emitted or absorbed by various plasma elements 
along the line-of-sight and having different oscillatory properties results in a noisy signal. 
After this leading work the subject of prominence oscillations remained dormant for about 
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15 years and revived in the 1980s, when a great amount of observations about prominence 
oscillations started being collected. 

According to the oscillatory amplitude, these events have been classified into two groups 
( [Oliver and Ballester|200"2"l l: 

- Large amplitude oscillations, that arise when the whole prominence is shaken by a 
Moreton wave (Moreton 1960) impinging on it (Ramsey and Smith 1966). As a con- 
sequence, the prominence gas undergoes a large displacement from its equilibrium po- 
sition and the prominence as a whole vibrates with considerable velocity amplitude (of 
the order of 20 km s or higher). These phenomena are reviewed in a separate work of 
this volume (Tripathi et al. 2009). 

- Small amplitude oscillations, with velocity amplitudes much smaller than those of large 
amplitude oscillations have also been frequently observed. The detected peak velocity 
ranges from the noise level (down to 0.1 km s _1 in some cases) to 2-3 km s -1 , although 
larger values have also been reported (Bashkirtsev and Mashnich 1984; Molowny-Horas 
|etal.|1999> . 

A distinguishing property of small amplitude oscillations, which gives them a distinction 
from large amplitude ones, is that they seldom influence the entire prominence. Thus, when 
spectrographic observations are carried out with a slit it is usually found that only a few 
consecutive points along the slit present time variations with a definite period, while all 
other points lack any kind of periodicity (e.g. Tsubaki and Takeuchi 1986; Balthasar and 



Wiehr 1994; Molo wny-Horas et al.||1997[ ). |Terradas et al.| ( |2002| l made a comprehensive 
study of the two-dimensional spatial distribution of Doppler velocity oscillations in a very 
large polar crown prominence and showed that, in their particular observation, oscillations 
are concentrated in a restricted area (54,000 x 40,000 km in size). Both propagating waves 
and standing oscillations are detected in this region and their wavelength and phase speed are 
determined. Using Hinode SOT observations of a quiescent prominence, Be rger et"aL| |2008 1 
also report on oscillations that do not affect the whole prominence body. These oscillations 
have periods from 20 to 40 minutes and typically last for only one or two cycles. They 
correspond to waves propagating vertically along threads with a projected phase speed of 
~ lOkms -1 . 

It has been mentioned that solar prominences are constituted by many thin, parallel 
magnetic threads filled with cold plasma, and as a consequence the dynamics of these com- 
ponents can form the basis of the dynamics of small amplitude oscillations. Early works fYij 
|et al.|199"T[|Yi a nd Engv oId|1991| l already noted the possible link between small amplitude 
prominence oscillations and the fibril structure. Unfortunately, the spatial resolution of the 
data analysed by Terrad aTet al.| ( |2002| > is not good enough to distinguish the prominence 
threads. It was necessary, hence, to wait until the advent of telescopes with much better spa- 
tial resolution to have observations in which the prominence fine structure is well resolved 
l |Lin et al.|2003|[Lin|2004[l. In th e analysis of the Doppler velocity in two threads belonging 
to the same filament, |Lin| ( |2004"} finds a clear sign of propagating waves and determines their 
period, wavelength and phase speed. This study is followed by a much more profound one 
in which the two-dimensional motions and Doppler shifts of 328 features (or "blobs") of 
different threads are examined. These features are observed to flow along the filament axis 
while oscillating at the same time. To simplify the examination of oscillations, |Lin| (2004 1 
computed average Doppler signals for each fibril and found that groups of adjacent threads 
oscillate in phase with the same period. This has two consequences: first, since the period- 
icity is outstanding in the averaged signal for each thread, the wavelength of oscillations 
is larger than the length of the thread. Second, fibrils have a tendency to vibrate bodily, in 
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groups, rather than independently. Ha observations conducted with the Swedish 1-m Solar 
Telescope by |Lin et alT] < |2007[ ) lead to similar results concerning the collective dynamics 
of fibrils, although propagating Doppler velocity signals with various periods and wave- 
lengths in other threads of the same filament are also detected. All these observations seem 
to indicate that prominence fibrils sometimes support collective oscillations and sometimes 
oscillate on their own. This topic deserves a more detailed observational study and, given 
the simplicity of fibrils compared to the full filament structure, theoretical investigations can 
give rise to fruitful comparisons with observations. 

During the 1980s some observational works, e.g. those by Tsubaki and Takeuchi ( 1986); 
Wieh r et aL] ( [T9 89 ), provided evidence about a temporal decrease of the oscillatory ampli- 
tude, which suggests that prominence oscillations are attenuated in time. A very clear exam- 
ple of this phenomenon can also be found in Figure 5a of |Landman et al.| ( |1977[ ). Here the 
integrated line intensity of the He D3 line displays oscillations with a period around 25 min 
and with an amplitude that decreases in time. These oscillations are present for a few cycles 
but, unfortunately, a precise value of the attenuation rate is not computed in this case (Fig- 
ure 3 of Berger et al. 2008 , also presents similar results, although with a much better spatial 
and temporal resolution). More detailed observational analyses of damped prominence os- 
cillations can be found in Molowny -Horas et al.| l |1999"| >; |Terradas et al.| ( |2002| l, who studied 
the spatial and temporal features of Doppler velocity oscillations in limb prominences. In 
the second of these works, already mentioned before, a sinusoidal function multiplied by 
a factor exp(— tj%) is fitted to the Doppler series, which allowed these authors to obtain 
values of the damping time, T, which are usually between 1 and 3 times the corresponding 
oscillatory period. In spite of the lack of similar studies, the existing evidence suggests that 
small-amplitude oscillations in prominences are excited locally and are damped in a few 
periods by an unknown mechanism. 

Many previous observational and theoretical results on prominence oscillations have 
been re viewed by e.g.|01iver and Ballester]l[27X)2^ ; |Engl^ ( [2l^ 

l [2006l l; [Banerjee et al.| ( |2007| >; |Engvold|p008| l, so our purpose here is to extend these review 
works by summarising other papers not discussed by these authors. Thus, we describe the 
first seismological analysis of prominence oscillations (Section [2| and the theoretical stud- 
ies on the attenuation of small amplitude oscillations (Sections[3][4]and[5]l. The plausible 
theoretical explanations of this effect are reviewed here and the relevance of different mech- 
anisms is assessed on the basis of the ratio of the damping time to the period (t/P), that 
must be in the range 1-10 to agree with the observational results. Works that invoke wave 
leakage from the prominence as the cause for the observed damping are described in Sec- 
tion [3] attenuation caused by non-ideal effects is discussed in Section [4] and wave damping 
by resonant absorption is presented in Section[5] Finally, Section [6] contains the concluding 
remarks. 



2 Seismology of active region threads observed with Hinode 

Observations with the Hinode Solar Optical Telescope by Okamoto et al. (2007) display an 
active region prominence composed of many fine horizontal threads that oscillate indepen- 
dently in the vertical direction as they move along a path parallel to the photosphere. Six 
different threads are studied and in each of them several points are considered. All points in 
each thread are found to oscillate in phase. The relevant parameters for a seismological anal- 
ysis of these events are flow velocities in the range 15^4-6 km s , oscillatory periods in the 
range 135-250 s and thread lengths in the range 1700-16,000 km. A prominence thread is 
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Fig. 2 Seismology of active region threads observed with Hinode. Thread (\>a p ) versus coronal (va c ) Alfven 
speed in two of the threads studied by Okamoto et al. ( 2007 ). The magnetic tube length ranges from 100,000 
km (bottom curve) to 250,000 km (top curve) in steps of 50,000 km. Along each curve the density ratio of 
thread to coronal plasma is 5 (*), 50 (0), 100 (A), 200 (□). Figures taken from |Terradas et al.|(2008) . 



a cold plasma condensation that occupies a segment of a much longer magnetic tube. In the 
present case it is not easy to directly measure the length of such magnetic tubes, although 
Oka moto et al.| ( |2007| l estimate that the wavelength of the oscillations is at least 250,000 
km and so the minimum length of the magnetic tubes is 125,000 km. |Terradas et al.1 (2008 1 
provide another estimation of this quantity, which they argue must be at least 100,000 km. 

These authors perform an independent seismological study of the six events by assuming 
that the thread is a dense plasma moving along a horizontal and straight magnetic tube which 
is tied to the dense photosphere at its ends. Additionally the low-jS and linear approximations 
are also imposed. For each thread three different studies are done; in increasing order of 
complexity they are: (a) thin tube approximation and no flow; (b) thin tube approximation 
and flow; and (c) full ideal MHD equations with flow. While the thin tube approximation is 
well supported by the thinness of the threads compared to their length, ignoring the effect 
of the flow does not seem a good idea because the detected flow velocities are rather large. 
Nevertheless, as we describe below the results prove just the contrary. 

In the absence of flow, case (a), transverse motions of the kind observed by Okamoto 
et al. (2007) can only be produced by the kink mode (Dymova and Ruderman 2005). The 
dispersion relation of the kink mode in a magnetic tube containing prominence material can 
be deduced from Equation (27) of Dymova and Ruderman (2005). Terradas et al. (2008) con- 
clude that, if no assumptions about the magnetic field, density and magnetic tube length (L) 
are made, then this dispersion relation contains too many unknowns and cannot be uniquely 
solved. But if a value of L is assumed (with the restriction L > 100,000 km from obser- 
vations), then a dependence of the prominence versus the coronal Alfven speed can be ob- 
tained. Figure [2] displays such a dependence for two of the six threads and several magnetic 
tube lengths. This figure shows that, unless the ratio of prominence to coronal density is 
unrealistically low (e.g. below 50), the Alfven speed in the thread reaches a minimum value 
that is not modified by that density ratio. This lower limit of va p is around 120 km s for 
thread #6 and lies between 200 and 350 km s _1 for the other five threads. 

To test the validity of this seismological result, Terradas et al. (2008 ) consider case (b), 
which involves solving a simple partial differential equation. Their numerical simulations of 
a thread that moves bodily along the magnetic tube yield in-phase oscillations of the whole 
thread, so a kink-like behaviour is recovered even in the presence of flow. These results are 
in excellent agreement with the observations of Okamoto et al. ( 2007 1. Perhaps surprisingly, 
for all six threads the period is only slightly smaller than that of case (a), the difference 
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being below 5%, i.e. smaller than the error bars of observations (of the order of or larger 
than 6.8%). Hence, the period of these transverse oscillations is almost insensitive to the 
thread flowing motions. Finally, Terradas et al. (2008) also investigate case (c) and recover 
the results of case (b), which implies that the thin tube approximation is quite reasonable, as 
expected. After these two studies the lower bounds for the Alfven speed derived from case 
(a) are firmly established and although this investigation does not lead to particular values 
of physical parameters, such Alfven speeds agree with the intense magnetic fields and large 
densities usually found in active region prominences. 



3 Damping of prominence oscillations by wave leakage 

Wave leakage is a common phenomenon by which the energy imparted by a disturbance 
to a given structure (e.g. a coronal loop or a solar prominence) is not confined as standing 
oscillations of the structure, which emits waves carrying energy into the external medium. 
This situation has been studied at length in relation with coronal loop oscillations; see Cally 
( 1986, 2003 ) and references therein. Most of the works described in this section make no 
mention of wave leakage, but our interpretation is that the obtained attenuation of oscilla- 
tions is simply a consequence of the emission of waves from the prominence into the corona. 



3.1 Line-current filament models 

In the celebrated model of Kuperus and Raadu|( [T974| ) a prominence is treated as a horizontal 
line current suspended in the corona. Such an approximation is justified by the fact that if 
the prominence plasma is supported against gravity by the Lorentz force, then a current runs 
along the filament. In this work the photosphere is substituted by a rigid, perfectly conduct- 
ing plane with a surface current distribution caused by the coronal magnetic field generated 
by the filament current. Such induced currents in turn give rise to a coronal magnetic field 
whose effect on the prominence is to provide the lifting force necessary to counteract grav- 
ity. The full magnetic field arrangement in the corona is determined by both the line current 
and the photospheric currents and can be of two types: either of normal polarity (NP) or of 
inverse polarity (IP); see Schutgens ( 1997b|; [van den Oord et al.| ( |199 8) for more details. 

This magnetic configuration was used by |van den Oord and Kuperus| ( |1992| >; |Schutgens| 
( |1997a|b[ l; |van den Oord et aT] ( |1998| l to study the oscillations of a filament. It is worth 
to remark that in these works the prominence is treated as an infinitely thin and long line, 
so that it has no internal structure, although the interplay of the filament current with the 
surrounding magnetic arcade and photosphere is taken into account. In addition, both NP 
and IP structures were considered; see Figure [3] in which the poloidal magnetic field is 
displayed (the toroidal magnetic field component being zero). The fundamental ingredient 
of these papers is that if a disturbance causes a displacement of the whole line current, 
that remains parallel to the photosphere during its motion, then the coronal magnetic field 
is reshaped and the photospheric surface current is modified. This, in turn, influences the 
magnetic force acting on the filament current. Such a force can either strengthen the initial 
perturbation, so that the original equilibrium is unstable, or diminish it, so that the system 
is stable against the initial disturbance. This simple description becomes more complicated 
when one takes into account that perturbations travel at a finite speed (namely the Alfven 
velocity) which causes time delays in the communication between the line current and the 
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Fig. 3 Magnetic field configuration of (a) normal polarity and (b) inverse polarity line-current prominence 
models used to investigate the damping of oscillations by wave leakage. The prominence electric current lies 
on the dip of magnetic field lines at the centre of the arcade in (a) and at the center of the magnetic island 
in (b). The magnetic field component perpendicular to the plane of the figure is zero. Figures taken from 
|Schutgens|jl997b) . 
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Fig. 4 Attenuation of prominence oscillations by wave leakage, (a) Quality factor (go = itt/P) of stable IP 
(solid curves) and NP (dashed curves) prominence oscillations as a function of the coronal Alfven speed, (b) 
Quality factor of the horizontally (squares) and vertically (diamonds) polarised stable oscillations versus the 
coronal density. Figures taken from Schutgens ( 1997b); Schutgens and Tdth 1 1999 1. 



photosphere. | van den Oo rd and Kuperus] ( | 1 992\ ; [Schutgens ( 199 7a|b[ l; |van den Oord et al.| 
( |1998| ) investigated the effect of these time delays on the filament dynamics. 

Exponentially growing or decaying solutions were found both for NP and IP prominence 
models; see Figure 4 of van den Oord and Kuperus ( 1992) and Figure 5 of Schutgens ( 1997a) 
as examples. An important conclusion that can be extracted here is that the attenuation is 
very efficient for many parameter values, and consequently the ratio of the damping time 
to the period is between 1 and 10 (i.e. in agreement with observations). This point is well 
illustrated by Figure [4^, in which the quality factor is depicted as a function of the Alfven 
speed. Different curves in this diagram correspond to different sets of parameters, with solid 
and dashed lines used to distinguish between IP and NP configurations. All the IP curves 
and one of the NP curves in this figure take values below Qq = 30, equivalent to t/P ~ 10, 
which is an indication of the efficiency of the damping of oscillations. 



3.2 Finite-thickness filament model 



An IP magnetic configuration similar to that used by |van den Oord a nd Kuperus ( 1992); 
|Schutgens||l997a|bfr ; |van den Oord et al.| ( |1998| ) was taken by Schutg ens and Toth 1 1999 i, 
although in this study the prominence is not infinitely thin and instead is represented by 
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a current-carrying cylinder. Schutgens and Toth (1999) carried out numerical simulations 
of the isothermal magnetohydrodynamic (MHD) equations and this implies that the equi- 
librium temperature is set to a constant value (10 6 K) everywhere. The authors mention 
that, despite the large temperature difference between solar prominences and the corona, 
this isothermal assumption can only have a minor effect. In addition, the photosphere is 
described as a perfectly conducting boundary, as in the papers examined above. The inner 
part of the filament is disturbed by a perturbation with a velocity amplitude of 10 km s 
and at an angle of 45° to the photosphere. This causes the prominence to move like a rigid 
body in the corona, both vertically and horizontally, and to undergo exponentially damped 
oscillations; see Figure 7 in Schutgens and Toth (1999). It is found that the horizontal and 
vertical motions are decoupled from one another (and so can be investigated separately) 
and that the period and damping time of horizontal oscillations are much larger than their 
respective counterparts for vertical oscillations. Again, strong damping can be achieved for 
some parameter values; see Figure [4]). It turns out that vertical oscillations are very effi- 
ciently attenuated for all the parameters considered in this work and that the same happens 
with horizontal oscillations for coronal densities above ~ 5 x 10~ 13 kg m~ 3 . These con- 
trasting properties of damped horizontal and vertical oscillations contain some potential for 
performing seismology of prominences. 



Now we turn our attention to the interpretation of the results presented so far. None of the 
works cited in the first part of this section, i.e. |van den Oord a nd Kuperus (1992); Schutgens 
( |1997a|b| l; |van den Oord et aT] ( |1998[ l, mention wave leakage as the cause of the attenuation, 
but instead these papers connect this behaviour with the effect of the distant photosphere 
and the presence of time delays in the communication of disturbances between the filament 
and the photosphere. On the other hand, Schutgens an d T6th| ([1999 ) link the damping of 
oscillations to the emission of waves by the prominence: the damping of horizontal motions 
is attributed to the emission of slow waves, whereas fast waves are invoked as the cause of 
the damping of vertical motions. Given that the main difference between this work and the 
other ones lies essentially in the cross section of the filament, there is no reason to believe 
that the physics involved are much different regardless of the prominence being modeled 
either as a straight and infinitely thin current or as a current-carrying cylinder. This issue 
must be examined in more detail to better understand the dynamics of solar prominences. 



A final remark about the exact nature of the wave leakage found by Schutg ens and T6th| 
( |1999| must be made. In this work (see their Figure 6a) the plasma-/3 in the prominence 
ranges from /3 > 1 in its central part to jS < 0. 1 at its boundary. Hence, waves emitted by the 
prominence into the corona propagate in a f5 <C 1 environment in which magnetic field lines 
are closed. Under these conditions, slow modes propagate along magnetic field lines and 
are unable to transfer energy from the prominence into the corona and so wave leakage in 
the system studied by Schutgens and Toth ( 1999) is only possible by fast waves. Then, it is 
hard to understand how the prominence oscillations can be damped by the emission of slow 
waves in this particular model, in which the dense, cool plasma is only allowed to emit fast 
waves. It must be mentioned, however, that the plasma-/3 in the corona increases with the 
distance from the filament, which implies that the emitted fast waves can transform into slow 



waves when they traverse the /3 ~ 1 region. This effect has been explored by McLaughlin 



|and Hood| ( |2006l l; |McDougall and Ho"od| l |2007l l, but see also references therein for similar 
studies. 
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4 Damping of small amplitude prominence oscillations by non-ideal effects 

Ballai (2003) took into account some dissipative mechanisms and explored, through some 
order-of-magnitude calculations, their importance as possible damping agents. In this work 
we are warned about the simplicity of the calculations and about the neglect of the solar 
corona, whose role on the damping of prominence oscillations is beyond this simple study. 
Dissipative mechanisms are separated into isotropic and anisotropic. Amongst the first ones, 
viscosity and magnetic diffusion are found to be rather inefficient in attenuating perturba- 
tions in only a few periods, as demanded by observations. Moreover, radiative losses (here 
modeled by means of Newton's cooling law, that gives an extremely simplified account of 
radiation in a solar prominence), may be of some importance, although their effect is diffi- 
cult to establish unambiguously because of the presence of the unknown radiative relaxation 
time-scale in Newton's cooling law. Therefore, radiation by the prominence plasma requires 
a more specific treatment to determine its significance. Regarding anisotropic mechanisms, 
in Ballai (2003 ) it is determined that, because of the very large density and very low temper- 
ature of prominences, viscosity can be considered isotropic and thus irrelevant as the cause 
of the observed damping of oscillations. In addition, thermal conduction is dominated by 
the presence of a magnetic field and is essentially parallel to the field direction. This mech- 
anism is very efficient as a damping agent for very short wavelengths. Since Ballai's work 
has been extended in other papers discussed below, the meaning of "very short wavelengths" 
will become clear later. Although other damping mechanisms are not contemplated in |Ballai| 
(2003}, e.g. wave leakage, ion-neutral collisions or resonant absorption, this work is impor- 
tant since it allowed subsequent investigations to discard the irrelevant mechanisms and to 
concentrate only in the rest. 



4. 1 Radiative losses and thermal conduction in a uniform medium 

The efficiency of thermal conduction (parallel to the magnetic field) and radiation in trans- 
porting heat, and thus in causing the attenuation of oscillations, has been addressed in a 
number of p apers; see|Ibanez S. and Escalona T] ( |1993) >; |Carbonell et al.| ( [2004 \ ; |Terr adas 



|et al.| ( |2003] l; |Carbonell et al. 1 2006 1. These works are concerned with the spatial and tem- 
poral damping of fast and slow waves propagating in a uniform medium. Radiative losses 
are simulated by Equation (7) of |Hildner| l |1974| l, which contains two parameters that can 
be tuned to represent optically thin or thick radiation; see |Carbonell et al.| ( |2004| l. Differ- 
ent sources of plasma heating are also included, although this mechanism has a negligible 
consequence on the properties of the waves. 

To understand the results of these studies it is worth to consider the characteristic time- 
scales of thermal conduction and radiative losses (t c and T r ), that following e.g. |De Moortel| 
|and H ood ( 2004 1 can be defined as 

- - LP CD 



(y-i)«i7" 



*r= 77. "„ 2 ^„ - (2) 



TP 

(y-l)p 2 XT a 

In these expressions p, p and T are the plasma pressure, density and temperature, L is 
a characteristic length-scale (the wavelength of oscillations, for example), Mi the parallel 
thermal conductivity and % and a are the two parameters in Hildner's radiative loss function. 
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Fig. 5 Time scales and critical wavelength in a 
uniform plasma with thermal conduction and ra- 
diative losses. Solid line: conduction time-scale, 
T c ; dashed line: radiative time-scale, l r \ dash- 
dotted lines: critical wavenumber that separates 
adiabatic from isothermal perturbations. The pa- 
rameter values used correspond to a prominence 
plasma: p = 5 x lO" 11 kg m~ 3 , T = 8000 K, 
fl = 0.8 (mean atomic weight), B = 10 G. More- 
over, parallel propagation (cos = 1) has also been 
assumed, so that for other propagation angles the 
critical wavenumbers are displaced to the right 
with respect to the vertical dash-dotted lines shown 
here. 



We see that T r is independent of L and so the rate at which the plasma looses energy through 
radiation does not vary with the wavelength or frequency of the waves. On the other hand, 
the conduction time-scale, T c , decreases linearly with L and so conduction becomes more 
effective at transporting heat for short wavelengths, just as noticed by Ballai (2003 ), for 
example. In the limit of very short wavelengths, conduction is so efficient that the isothermal 
regime is reached. The critical wavenumber, k c , at which this transition occurs is given by 
|Porter et al.| ( fT994| > as 

_ 2nk B c s 2nk B v A 

"■c—slow — i ^c— fast — 9 - ; \P) 

Mi cos Q Mil cos^ 8 

for slow and fast waves, respectively. Here n is the number density, k B Boltzmann's constant, 
c s and \'a the sound and Alfven speeds and Q the angle between the wavenumber and the 
unperturbed magnetic field. To derive these formulas it is necessary to assume c s <C Va, 
together with (Or w k x c s for the slow mode and (Or rs kv A for the fast mode, where (Or 
is the real part of the frequency and k x the component of the wavenumber parallel to the 
equilibrium magnetic field. 

Figure[5]presents a plot of the characteristic conduction and radiative times together with 
the critical slow and fast mode wavenumbers marking the transition from the adiabatic to the 
isothermal regime. Now, let us consider a perturbation with tunable wavelength travelling 
in a magnetised, uniform medium in which thermal conduction and radiation can transport 
heat. Let us start with an extremely short wavelength, so that the wavenumber is larger than 
k c and the perturbation is isothermal (hence, we are at the right of Figure |5j. Let us now 
increase the wavelength until the wavenumber becomes just smaller than k c . Now perturba- 
tions are no longer isothermal, but given that thermal conduction is still more efficient than 
radiation, since it works in a very short time-scale (i.e. T c <gC T r ), the perturbation is domi- 
nated by conduction. Next the wavelength is allowed to increase further and consequently 
the conduction time-scale increases proportionally to L. Under these conditions, conduction 
becomes less and less efficient and there is eventually a wavenumber (or length-scale) at 
which the conduction and radiation time-scales become equal; it corresponds to the position 
in Figure [5] in which the solid and dashed lines cross. This length-scale, L*, is obtained by 
imposing T c = T,- and with the help of Equations (|T]| and <j2j) is 



(4) 



p 2 xT a-l- 

Finally, for wavelengths larger than L* radiative losses dominate the energy budget of the 
plasma. 
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Fig. 6 (a) Damping time and (b) damping per period, D p = 2%% /P, versus the wavenumber of the fast 
mode in an unbounded prominence medium with thermal conduction and radiative losses. The three curves 
correspond to different parametrisations of the prominence radiative losses. The wavenumber is given in m~ 1 . 
Figures taken from Carbonell et al. (2 004) . 

These formulas are useful to explain some features of the damping time (in the case 
of temporal damping of oscillations) or the damping length (in the case of spatial damping) 
obtained by|Ibanez S. and Escalona T.|fl993l l; |C"arbonell et al^ l [2T)04) ; |Terradas et al.| (|2005 ); 
|Carbonell et al.| ( |2006| l. Consider, for example, the damping time and the damping per period 
{D p = 2ni /P) of the fast mode modified by radiative losses and thermal conduction in a uni- 
form medium. Their variation with the wavenumber is displayed in Figure|6] in which each 
of the three curves originates from a particular pair of values of % and a, the two parameters 
in the radiative loss function. In our discussion we concentrate on the solid line, although 
similar conclusions can be deduced for the other two (and also for the slow mode, which is 
not examined here). All curves in Figure [6] display three changes of slope whose origin can 
be partly explained on the basis of the transitions between different regimes studied above. 
For the paramater values used in this work, we have k c -f ast = 1.82 m and L* = 4770 m, 
that corresponds to the wavenumber £:* =2k/L* = 1.32x 10~ 3 m~ 1 . The right-most change 
of slope in Figure[6]is around logk = 0.5 and coincides roughly with the value of k c -f ast , so 
it is caused by the transition between the isothermal and the conduction-dominated regimes. 
In addition, the intermediate change of slope is around log A: = — 3 and its agreement with k* 
confirms that it corresponds to the transition from the conduction-dominated to the radiation- 
dominated regimes. Therefore, the shift from one regime to another leaves its imprint in the 
variation of T and x/P with respect to the wavenumber. We also note that, since the dif- 
ference between the three curves of Figure [6] lies in the radiative loss function used, they 
overlap in the conduction dominated regime. 

The main conclusions of Carbonell et al. (2004); Terradas et al. (2005); Carbonell et al. 
( 2006 ) regarding the damping of prominence oscillations are that the slow mode is strongly 
attenuated for wavenumbers within the range of observed values (i.e. for k in the range 
10~ s — 10~ 6 m ), while the fast mode is not so strongly damped and its attenuation rate is 
too small compared with the values from observations. Moreover, the considered approxi- 
mations for optically thin or thick plasmas can give very different attenuation properties and 
all the considered heating mechanisms are of no relevance. 



4.2 Radiative losses and thermal conduction in a structured medium: prominence slab 

The influence of the surrounding corona on the damping of magnetoacoustic oscillations via 
thermal conduction and radiative losses has been studied by |Soler et al.| (2007a 2 009a| ). 
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Fig. 7 Sketch of a prominence slab (temperature T„ and density p p ) embedded in the solar corona (tempera- 
ture T c and density p c ). Left: longitudinal magnetic field, parallel to the prominence length. Right: transverse 
magnetic field, perpendicular to the prominence length. Figures taken from Soler et al. 1 2007a. 2009al. 




Fig. 8 Phase speed (in km s~') versus the di- 
mensionless longitudinal wavenumber of the nor- 
mal modes of a prominence-corona system with 
longitudinal magnetic field. Solid lines denote 
kink modes and dotted lines denote sausage 
modes. Only the fundamental modes are plotted. 
The shaded zones are forbidden regions where 
evanescent-like solutions in the corona do not ex- 
ist. Note that the vertical axis is not drawn to scale. 



Figure taken from Soler et al. ( 2007b) . 



In |Soler et ah] p007a| l the filament is taken as a plasma slab embedded in the solar 
corona and with a uniform magnetic field parallel to the filament axis; see Figure [7^. Both 
the prominence plasma and the coronal plasma are uniform and threaded by the same mag- 
netic field. To understand the results of Sol er et aE] ( |2007a| l it is necessary to describe the 
main properties of the linear, adiabatic normal modes of such a configuration, that have been 
studied by |Edwin and Roberts| ( fl982) ; |Joarder and Roberts] ( |1992a| ; [Soler et al.| l (2007b^ In 
the latter work the solutions are classified in three types (Figure [8b: fast, internal slow and 
external slow modes and it is concluded that the perturbations associated to internal slow 
modes are essentially confined to the prominence slab, so that these waves are hardly influ- 
enced by the coronal environment. On the other hand, fast modes have tails that penetrate 
in the corona and as a result coronal conditions are important in determining the features of 
these waves. Moreover, the confinement of fast modes becomes poorer for small values of 
the wavenumber in the direction of the filament axis. Regarding external slow modes, they 
produce almost negligible perturbed amplitudes inside the cold slab and thus seem uninter- 
esting in the context of prominence oscillations. 

Soler et al. ( 2007a) included the thermal conduction, radiative losses and plasma heating 
terms in the energy equation and studied the damping properties of the fast and internal and 
external slow modes. They obtained values of t/P in agreement with observations both 
for the internal slow and the fast modes; this last result implies that the presence of the 
surrounding corona modifies the damping properties of fast waves. Such as expected, the 
damping rate of the internal slow mode is not modified by the inclusion of the corona in 
the model and it only depends on the value of the wavenumber in the direction parallel to 
the magnetic field. The attenuation of this mode can be recovered by inserting the parallel 
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wavenumber, k x , in the expression CO = k x c s and by substituting c, by a "modified sound 
speed", A, whose square is given by 

, c 2 (y~ (iH^ + tOr-Vp) + i 7 C0 

A 2 = — ^— -J- , (5) 

7 (r- 1 ) (]r«qi** +«»*)+»«> 

where co p and C0j are determined from the partial derivatives of the heat-loss function with 
respect to density and temperature, respectively. The damping time of the fast kink mode, 
however, is severely influenced by the addition of the coronal medium. This can be clearly 
appreciated by comparing the damping time in a uniform medium (Figure |6^) and in a 
structured medium with longitudinal magnetic field (Figure[9]l. The reason for such a strong 
modification of the fast kink mode damping time is, firstly, that this mode couples with the 
external slow mode, whose properties are essentially dictated by the coronal plasma. And 
secondly, apart from this coupling the coronal medium has a direct influence because of the 
presence of coronal tails in the perturbed variables. In Sol er et aT] ((2007a 1 a close exam- 
ination of the role of the corona in the damping of the fast mode was done by removing 
one of the important non-adiabatic mechanisms (thermal conduction or radiation) either in 
the filament or in the corona and then computing the period and damping time. The rele- 
vance of each mechanism is revealed by a discrepancy between the damping time computed 
with and without this mechanism (Figure |5J. The conclusion of this study is that coronal 
mechanisms govern the attenuation of the fast modes for wavenumbers k x < 10~ 7 m _1 , 
whereas prominence mechanisms prevail for k x > 10 m . In each of these ranges radi- 
ation is dominant in the lcw-k x range and thermal conduction is dominant in the high-^ 
range, such as found in Section |4~T| for a uniform medium. More specifically, coronal radia- 
tion controls the attenuation for k x < 10~ 9 m _1 , coronal conduction does so for 10~ 9 m _1 
< k x < 10~ 7 m _1 , while prominence radiation and conduction are the most important non- 
adiabatic mechanisms for 10~ 7 m _1 < k x < 10~ 3 m _1 and k x > 10~ 3 m _1 , respectively. 
Hence, in the range of observed wavelengths both thermal conduction in the corona and 
radiation from the prominence gas are important to explain the reported damping times. In 
Section |4~T1 the characteristic time-scales of thermal conduction and radiation in a uniform 
medium were invoked to explain the transition from the radiation-dominated wavenumber 
range to the conduction-dominated wavenumber range. The same idea can be applied to 
explain this transition in the present, structured configuration both for the prominence and 
for the corona. In the case of the prominence, the transition value is the same as before, 
namely k* = 1.32 x 10 -3 m _1 , whereas for the corona it is k* = 3.42 x 10~ 9 m . These 
numbers coincide well with the transitions at which different mechanisms dominate in Fig- 
ure]^ Moreover, the wavenumber at which the transition between the ranges of dominance 
of coronal conduction and prominence radiation takes place can be derived from the equality 
of T c in the corona and x r in the prominence. This calculation yields k = 1.32 x 10~ 7 m _1 , 
also in good agreement with the value from Figure [9] More details can be found in Soler 

|eTal~| i |2007a|l. 

In |SoIe7et1iI~| l |2009a| l 

a prominence-corona configuration with an equilibrium magnetic 
field perpendicular to the prominence slab is considered; see FigurejT}}. Such a structure can 
support fast and slow modes of internal or external character, depending on whether their 
properties are governed mainly by the prominence or the corona, respectively. Moreover, 
there is another type of fast and slow mode, called string by Joarder and Roberts ( 1992b I 
and hybrid by jOliver et al.|(T993) , that is simultaneously internal and external. |Soler et al.| 
(2009a ) find that slow modes are efficiently damped by the considered non-adiabatic mech- 
anisms (thermal conduction and radiative losses), such as found by |Soler et al.| ( |2007a[ l for 
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Fig. 9 Damping time versus the longitudinal 
wavenumber for the fundamental fast kink mode 
in a prominence slab with longitudinal magnetic 
field. Different linestyles represent the results ob- 
tained when different non-adiabatic mechanisms 
are omitted: all mechanisms considered (solid 
line), prominence conduction eliminated (dotted 
line), prominence radiation eliminated (dashed 
line), coronal conduction eliminated (dot-dashed 
line) and coronal radiation eliminated (three dot- 
dashed line). The shaded region indicates the in- 
terval of observed wavelengths. Figure taken from 
|Soleretal.|j2007a} . 



a perpendicular orientation of the magnetic field. It is also shown that prominence radiation 
dominates the attenuation of internal slow modes and that prominence radiation and coro- 
nal conduction together dominate that of the hybrid slow mode. This result for the hybrid 
mode is coherent with the fact that its perturbations achieve large amplitudes both in the 
prominence and the corona, so one expects that the most relevant damping mechanisms of 
each medium govern together the attenuation of the hybrid mode. On the other hand, fast 
modes are very poorly damped in the present configuration, which implies that the assumed 
magnetic field orientation is decisive for these modes. A remarkable feature of fast modes 
is that they can become unstable under certain circumstances. Thermal conduction allows 
energy transfer between the prominence slab and the coronal medium. Prominence radia- 
tion has an essential role in dissipating the extra heat injected from the corona and stabilises 
the oscillations. Nevertheless, thermal instabilities appear if radiative losses from the promi- 
nence plasma are omitted or significantly reduced (e.g. caused by an increase of the optical 
thickness) since the plasma cannot dissipate the extra injected heat in an efficient way. Given 
that different parametrisations of the radiative loss function have been given in the literature 
( |Hildner|1974]|Rosner et al.|1978"}|Mime et al.|1979| > and that they correspond to different 
optical thicknesses, the relevance of the optical thickness in the stability of fast waves is not 
firmly established yet. 



4.3 Attenuation by radiative losses and thermal conduction in the presence of flow 

Material flows are common features in solar prominences ( |Zirker et al.|1998[|Lin et al.T2 003 
2007). Carbonell et al. (2009) investigated the attenuation of the slow and thermal waves by 
radiative losses and thermal conduction in a uniform medium with flow. In this work the fast 
mode is ignored because these non-adiabatic effects do not contribute in a significant way to 
its attenuation ( Carbonell et al. 2004 1; slow and thermal waves are investigated in the limit of 



field-aligned propagation. Carbonell et al.| ( [2009| ) find that the period and damping per period 



can show a strong dependence on Vo, the flow speed, for certain values of this parameter; 
Figure [TO] presents their results for the slow mode. They also find that the greatest period 
and damping per period of slow waves is obtained for flow speeds close to the real part of 
the non-adiabatic sound speed, which in the presence of flow is deduced from Equation |5]i 
with co substituted by CO — k x VQ. Regarding the thermal mode, which in the absence of flow 
does not propagate, it is found that it becomes propagating for vo 7^ 0. Therefore, when 
flows and time damped oscillations are observed in a prominence/filament, this wave should 
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be also taken into account as potentially responsible of the observed damped oscillations. 
These authors remark that it may not be possible to determine whether an observed period 
and damping time are associated to a slow wave or a thermal wave in a flowing material 
and that a possibility to determine the origin of such wave is to resort to the temperature 
perturbation, which should be rather large for the thermal wave and much smaller for the 
slow wave. 
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Fig. 10 (a) Period and (b) damping per period, D p = 2nz/P, versus the flow speed for slow waves in a 
uniform plasma with flowing motions and subject to radiative losses and thermal conduction. Solid, dotted, 
dashed and dash-dotted lines correspond to wavenumbers k x = 10~ 8 , 10~ 7 , 10~ 6 , 10~ 4 m~', respectively. The 
inclusion of flow in the model gives rise to the presence of two slow modes for each wavenumber when the 
flow speed is non-zero. Figures taken from Carbonell et al. 12009). 



4.4 Radiative losses and thermal conduction in a structured medium: cylindrical 
prominence thread 
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Fig. 11 (a) Sketch of the thread model used bylSoleretaL 2008 to study the damping of oscillations caused 
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by thermal conduction and radiative losses in the presence of mass flow, (b) Damping time of the fast kink 
mode in a cylindrical thread (solid) and in a thin slab with longitudinal magnetic field (dotted). In both cases 
the flow velocity is zero. The shaded area indicates the interval of observed wavelengths. Figures taken from 
|Soler et al.|(20(%V 



Instead of considering the attenuation of perturbations in a full-prominence configura- 
tion, as in Section |4~2] S oler et aI7| ([2008i) restrict themselves to a single thread. Their model 
consists of a homogeneous, infinite cylinder representing the thread embedded in an un- 
bounded corona; see Figure [TTfc . Two additional ingredients of the model are a uniform 
magnetic field both inside and outside the thread and flow motions parallel to the magnetic 
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field in both regions. Thermal conduction and radiative losses are taken into account as 
damping mechanisms, and the combined effect of these non-ideal physics and the steady 
flow on the attenuation of oscillations is assessed. In this work it is found that, in the ab- 
sence of flow, slow modes are efficiently damped by non-adiabatic effects, while fast kink 
waves are in practice not attenuated, since their damping times are much larger than typi- 
cal lifetimes of filament threads. Therefore, the damping by non-adiabatic mechanisms of 
kink oscillations is much less efficient in this cylindrical case than in the slab geometry 
(Figure [TTp). In the presence of flow it is well-known that the symmetry between waves 
propagating parallel or antiparallel to the flow is broken. The thermal mode behaves as a 
wave that propagates parallel to the flow, and its motions are mainly polarized along the 
longitudinal direction. Nevertheless, this oscillatory behaviour is not likely to be observed 
in practice due to its quick damping. The damping time of both slow and thermal waves is 
not affected by the flow. On the contrary, for realistic values of the flow velocity, the larger 
the flow, the larger the attenuation of parallel fast kink waves, whereas the contrary occurs 
for antiparallel fast kink solutions. Nevertheless, this effect is not enough to obtain realistic 
damping times in the case of fast kink modes. 



4.5 Partial ionisation effects in a uniform medium 

Chromospheric and prominence plasmas are not fully ionised and so several new effects are 
present in comparison to a fully ionised plasma. For example, electric charges are frozen 
to magnetic field lines, but neutrals are not and thus the neutral and ionised fractions of 
the plasma behave differently. Collisions between neutrals, on one hand, and electrons and 
ions, on the other hand, arise and consequently a modified Ohm's law is obtained. As a 
consequence, the main outcome of the interaction between neutrals and charged particles is 
the presence of enhanced Joule heating and enhanced magnetic diffusion. These ideas were 
presented by Piddington l |1956| l in connection with the heating of the solar chromosphere 
and corona and have also been invoked to explain the dynamics of chromospheric spicules 
(e.g. Haerende l|1992| James and Erdelyi|2002| > and of flux tubes emerging from the solar 



interior jArber et al.|2007| > 



The influence of partial ionisation of the prominence plasma on the properties of the 
fast and slow modes has been investigated by |Forteza et al.| ( |2007[ >. In this work the one- 
fluid MHD equations for a partially ionised plasma made of electrons, protons and neutral 
hydrogen are derived. Next, a uniform medium with a straight magnetic field is considered. 
The linear regime is assumed and this results in the vanishing of the Joule heating term in 
the energy equation, so that the corresponding physical effect is absent. One may conclude 
that in more realistic configurations dissipation may be stronger than that found by |Forteza| 
et al. (2007). The main result of this work is that the fast mode can be strongly damped 
(in agreement with observations) for small values of the ionisation fraction (Figure 12). 
whereas the slow mode is hardly attenuated. The mechanism responsible for the damping of 
the fast mode is Cowling's diffusion, that vanishes in a fully ionised plasma. An analytical 
approximation for the fast mode damping rate, 

h'-jr^ 1 *" (6) 



is derived by Forteza et al. ( 2007 1, with ^„ the relative density of neutrals, p the gas density, B 
the magnetic field strength and k the modulus of the wavenumber. Such a simple expression 
can be of practical use in the interpretation of observations. It indicates that the damping is 
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stronger, i.e. x/P is smaller, for almost neutral plasmas with low density, strong magnetic 
field and for short wavelengths. 
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Fig. 12 Damping rate versus the ionisation frac- 
tion, fl = 1/(2 — ^„), for the fast mode in a uni- 
form, partially ionised plasma. The solid line rep- 
resents the numerical solutions of the full disper- 
sion relation while the triangles represent the re- 
sults obtained by solving the aproximate Equa- 
tion J6j. Figure taken from Forteza et al. ( 2007 ). 



4.6 Radiative losses and thermal conduction in a uniform, partially ionised medium 
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Fig. 13 Ratio of the damping time to the period as a function of the wavenumber, k, for the (a) slow and (b) 
fast waves in a partially ionised plasma with thermal conduction and radiative losses. Different line styles 
correspond to different ionisation fractions: p. = 0.5 (dashed), fl = 0.8 (solid), jl = 0.95 (dash-dotted) and 
p. = 0.99 (dotted). The shaded region corresponds to the interval of observed wavelengths in prominences. 
Figures taken from Forteza et al. (20081. 



A natural extension of Carbonell et al. (2004 ), concerned with the attenuation of fast and 
slow waves by thermal conduction and radiation, and Forteza et al. (2007), concerned with 
the attenuation of the same waves by ion-neutral effects, is to consider together a partially 
ionised plasma with thermal conduction and radiative losses. This task was performed by 
|Forteza et al.| |2008), who derived the equations governing the dynamics of such medium. 
A particularly interesting issue is the different conduction properties of electrons (mainly 
along magnetic field lines) and neutrals (isotropic), which implies that the contribution of 
the neutral species to the thermal conduction must be taken into account separately following 
|Parker| ( [T953) ); |Ibanez and Mendoza| ( [19901 ). |Forteza et al.| < [2008| > find that the period of the 
magnetoacoustic waves remains basically the same as in the ideal case and that their model 
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gives values of the ratio of the damping time to the period similar to the ones obtained in 
observations for the three wave modes (fast, slow and Alfven). In addition, the inclusion of 
non-adiabatic terms in the partially ionised set of equations increases the damping of fast and 
slow waves in the interval of observed wavelegths as compared with the results obtained for 
a non-adiabatic fully ionised plasma ( Carbonell et al. 2004 1 and an adiabatic partially ionised 
plasma I Forteza et al.|2007 1. For slow waves, the minima of the ratio of t/P, corresponding 
to attenuation maxima, are displaced towards longer wavelengths as compared to when only 
thermal conduction and radiative losses are included in the model. An increase of the neutral 
portion in the plasma produces a displacement of these ranges of maximum damping to 
longer wavelengths (Figure [T3fr). Radiative losses are dominant for long wavelengths (just 
as found by Carbo nell et al.|20 04), while the rest of the wavenumber interval is dominated 
by thermal conduction and ion-neutral collisions. Regarding fast waves, radiation is the 
dominant damping mechanism for long wavelengths, while in the rest of the considered 
wavenumber interval the attenuation is dominated by the effect of ion-neutral collisions. 
An important result is that fast waves only exist for wavenumbers smaller than a critical 
value that depends on the ionisation fraction (see Figure [T3p). In spite of this, the critical 
wavenumber is large in comparison with the typical wavenumbers of waves in prominences. 
Finally, for typical prominence temperature values, the contribution of electrons to thermal 
conduction is negligible in front of the contribution of neutrals. 



5 Damping by resonant absorption in a prominence thread 




Fig. 14 Sketch of a radially non-uniform filament 
fine structure of mean radius a, with a transverse 
inhomogeneous layer of thickness /. For 1 = 
the geometry of the model reduces to that of Fig- 
ure | 1 1 fr . Figure taken from Arregui et al. 1 2008 1. 



The plasma density varies by about two orders of magnitude between a filament thread 
and the surrounding corona. In such a highly inhomogeneous configuration fast kink modes 
can be efficiently damped by transferring their energy to Alfven modes through resonant 
absorption. Arregui et al. (2008 ) examine the relevance of this mechanism in the attenuation 
of small amplitude prominence oscillations by modeling the magnetic and plasma configu- 
ration of an individual fibril as an infinitely long, straight, cylindrically symmetric flux tube 
of mean radius a surrounded by a coaxial layer of thickness / (Figure 14 1. The /3 = approx- 
imation is used. The coronal and thread densities are denoted by p c and pf, respectively, and 
their ratio as c = Pf/p c - This parameter is much larger than unity and probably of the order 
of 100 or larger. In this simple model, the magnetic tube representing the thread only con- 
tains cool material, i.e. the presence of hot material inside the tube is ignored. Taking into 
account that the observed propagating waves in threads have a wavelength, A, much larger 
than the tube radius, [Arregui et aL"| lf2008 ) start by using the long-wavelength approximation 
to make some analytical progress. Next, they use the so-called thin boundary approximation 
(/ <C a) to derive the following expression for the ratio of the damping time to the period 
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T -=F a - C -±± (7) 
P Ic-V {> 

where F is a numerical factor that depends on the particular variation of the density in the 



1988 



Goossens et al. 



nonuniform layer: for a linear variation F = A/n ( Hollweg and Yang 
1992), while for a sinusoidal variation F = 2/tt (Ruderman and Rob erts|2002| >. A simple 
substitution in this formula of a density ratio c = 200 and l/a = 0.1, which corresponds 
to a transitional layer whose thickness is 10% the thread radius, yields x/P ~ 6, in good 
agreement with observations. Then, this analytical approximation seems to indicate that the 
process of resonant absorption is an efficient damping mechanism for fast MHD oscilla- 
tions propagating in these structures. Next, Arregui et al. (2008) solve numerically the full 
set of linear, resistive, small-amplitude MHD wave equations for the transverse kink os- 
cillations. In this manner, the approximations that lead to Equation ^ are avoided. These 
numerical calculations confirm the analytical estimate that, for the typical large filament to 
coronal density contrast, the mechanism produces rapid damping on timescales of the or- 
der of a few oscillatory periods only. They also reveal that for thin layers (l/a ~ 0.1) the 
inaccuracy of the long- wavelength approximation produces differences of up to ~ 10% for 
the combination of short wavelength with high-contrast thread. For thick layers, differences 
of the order of 20% are obtained (in agreement with |Van T3o orsselaer e et al.|2004) . More- 
over, the authors conclude (as shown by Figure [T5jl that the obtained damping rates are only 
slightly dependent on the wavelength of perturbations and that the damping rate becomes 
independent of the density contrast for large values of this parameter. This last result has two 
seismological consequences. First, the observational determination of the density contrast is 
less critical than it is in the low-contrast regime, the one corresponding to coronal loops. 
Second, according to seismic inversion results that combine the theoretical and observed 
periods and damping times (Arregui et al. 2007; Gooss ens et al.|20 08), high-density thread 
models would be compatible with relatively short transverse inhomogeneity length scales. 
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Fig. 15 Ratio of the damping time to the period for fast kink waves in filament threads (radius =100 km) 
attenuated by resonant absorption, (a) x/P as a function of the density contrast with l/a = 0.2 for two longi- 
tudinal wavelengths, (b) x/P as a function of the transverse inhomogeneity length scale for two combinations 
of wavele ngth and density contrast. In both pan els solid lines corresp ond to the analytical solution given by 
Equation j!5| with F = 2/jl. Figures taken from Arregui et al. 12008). 



In a subsequent work, Soler et al. 1 2009b I note that in a j3 ^ plasma with prominence 
conditions embedded in the much hotter corona the fast kink mode can also undergo a reso- 
nant coupling to slow continuum modes. Hence, it is argued, the combined effect of the two 
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resonant couplings may lead to an enhanced attenuation rate of the kink mode. The struc- 
ture considered by |Soler et al.| ( |2009b| l is a thread identical to that of Figure[l4j although the 
jS = assumption is removed from their calculations and, consequently, the plasma pressure 
is finite. Such as in Arr egui et al.| ( |2008| l, in this work a combination of analytical expres- 
sions and numerical analysis is used. The analytical formulas not only allow to predict the 
positions of the slow and Alfven resonances, but they also provide the contribution of each 
of the two mechanisms to the damping time, Toa and Xd s , whose ratio is 



(8) 



with m and k z the azimuthal and longitudinal wavenumbers and c s and va the sound and 
Alfven speeds. Using m = 1 and a typical value for the wavelength (k z a = 10~ 2 ) results in 
Tda/tdj ~ 10 , which means that the resonant coupling of the kink mode with the Alfven 
continuum produces an attenuation rate that is 10 7 times smaller than that coming from the 
coupling with the slow continuum. These results are also verified from the numerical com- 
putations (see Figure [T6|, which confirm that resonant absorption with the slow continuum 
is irrelevant in the damping of the observed transverse thread oscillations. 
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Fig. 16 Ratio of the damping time to the period as 
a function of the dimensionless wavenumber, k z a, 
corresponding to the kink mode for / /a = 0.2. At- 
tenuation in this cylindrical thread model is caused 
by resonant coupling to Alfven and slow contin- 
uum modes. The solid line is the complete numer- 
ical solution, whereas the symbols and the dashed 
line are the results of the thin boundary approxi- 
mation corresponding to the Alfven and slow reso- 
nances, respectively. The shaded zone corresponds 
to the range of typically observed wavelengths in 
prominence oscillations. Figure taken from |Soler] 
|eTaLl(2009b) . 



6 Conclusions 

Although the study of small-amplitude waves and oscillations in prominences constitutes 
a discipline that might complement the direct determination of prominence parameters by 
providing independent values based in the comparison between observations and theory, 
so far prominence seismology is more a promise than a reality because of the large gap 
between observation and theory. The analysis presented in Section[2]is a good example of a 
seismological application to prominences, but, unfortunately, it is the only one of its kind. 
Such a gap arises because of the little restrictions imposed by observational works (which 
sometimes reduce to reporting the period of the detected oscillations) and the simplicity of 
theoretical studies (which neglect the intricate nature of prominences and substitute it by a 
very simplified physical model). Nevertheless, these two sides are coming together as the 
complexity of works increases. 
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In this paper the possible role of several physical mechanisms in the attenuation of 
prominence oscillations has been reviewed. It has been shown that both the fast and slow 
waves can undergo strong damping under a variety of physical conditions and with the in- 
tervention of different effects. Nevertheless, the list of the presumably relevant mechanisms 
has not been totally explored and other effects should be investigated. In addition, the geom- 
etry of the models is too simplified, which implies that the results derived so far only give a 
rough guide to the damping features in actual prominences. The efficiency of wave leakage 
and coronal conduction in draining oscillatory energy from the prominence and transferring 
it to the corona calls for the study of complex configurations in which the structure of the 
full prominence-corona system is treated more realistically. 

Most of the works are based on the linear approximation, whose validity seems quite 
robust in view of the small amplitude of prominence oscillations in many cases. Neverthe- 
less, apart from the so-called small-amplitude prominence oscillations, that in general affect 
small regions of a prominence and have Doppler velocity peaks typically below 1-2 km 
s _1 , there is a different kind of phenomenon, the large- amplitude oscillations (Tripathi et 
al. 2009, this volume). They are characterised by a much highe r amplitude (with oscillatory 
velocities up to 90 km 1 ) and disturb the whole filament; see Eto et al. |2002|; Jing et al. 



2003 



; |Okamoto efaL] < (2004) |; |Jing et al.| < [2006l >; |Isobe and Tripafhi| i |2006| >; |Vrsnak et al 



2007 1. These prominence vibrations also damp very rapidly, in a few periods, so it is in 
order to ascertain whether the mechanisms at work in the damping of small-amplitude os- 
cillations are also the main ones for the large-amplitude counterparts or whether non-linear 
effects can also cause the attenuation of the latter. 
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